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SENSITIVITY ANALYSIS FOR SHAPE OPTIMIZATION OF A 
FOCUSING ACOUSTIC LENS IN LITHOTRIPSY 

VANJA NIKOLIC AND BARBARA KALTENBACHER 


Abstract. We are interested in shape sensitivity analysis for an optimization problem 
arising in medical applications of high intensity focused ultrasound. The goal is to find 
the optimal shape of a focusing acoustic lens so that the desired acoustic pressure at a 
kidney stone is achieved. Coupling of the silicone acoustic lens and nonlinearly acoustic 
fluid region is modeled by the Westervelt equation with nonlinear strong damping and 
piecewise constant coefficients. We follow the variational approach to calculating the shape 
derivative of the cost functional which does not require computing the shape derivative of 
the state variable; however assumptions of certain spatial regularity of the primal and the 
adjoint state are needed to obtain the derivative, in particular for its strong form according 
to the Delfour-Hadamard-Zolesio Structure Theorem. 


1. Introduction 

The present paper is concerned with the shape optimization problem arising in lithotripsy, 
where an optimal focusing of the ultrasound waves is needed in order to concentrate the ul¬ 
trasound pressure at the kidney stone and avoid lesions of the surrounding tissue. 

The design of currently used devices in lithotripsy is mainly based on two different prin¬ 
ciples (see, e.g, [18], [19] , [20], [31]): excitation and self-focusing by a piezo mozaic or an 
electro-magneto-mechanical principle. 

Axis of Rotation 



based on the electromagnetic principle 

We investigate the latter case, where Lorentz forces acting on a membrane radiate an acous¬ 
tic pulse in a fluid; the pulse is then focused by a silicone lens at a kidney stone. The aim of 
obtaining a sharp focus of the acoustic pressure exactly at the desired location leads to the 
task of optimizing the shape of the acoustic lens. We will consider the case where the lens 
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is modeled as an acoustic medium surrounded by a nonlinearly acoustic fluid. 

One of the commonly used models for propagation of nonlinear ultrasound is the Wester- 
velt equation 

(1.1) (1 — 2 ku)u — c 2 Au — bAu = 2 k(u) 2 , 

given here in terms of the acoustic pressure u, where b denotes the diffusivity and c the speed 
of sound, k = /3 a /X, A = gc 2 is the bulk modulus, g is the mass density, /3 a = 1 + B/(2A), 
and B/A represents the parameter of nonlinearity. For a detailed derivation of (11.11) we refer 
to m, B29 and [32]. 

Westervelt’s equation is a quasilinear wave equation which can degenerate due to the 
factor 1 — 2 ku. To avoid this degeneracy, i.e. to find an essential bound for the acoustic 
pressure u, Sobolev embedding H 2 { f2) —>• L°°(Q) is employed (cf. [15], [IS]) and thus for 
this model only very smooth solutions can be shown to exist. However, H 2 regularity on the 
whole domain is too high of a demand in the case of interface coupling of acoustic regions 
with different material parameters. 

In [4], a nonlinear damping term was introduced to the Westervelt equation 

(1.2) (1 — 2 ku)u — c 2 Au — div(b((l — S) + (5|Vh|' ? ” 1 )Vu) = 2 k(u) 2 , 

with 6 € (0,1), q > 1, q > d — 1, d € {1,2,3}, which allowed to show existence of weak 
solutions with W l,q+l regularity in space, and in turn well-posedness of the acoustic-acoustic 
coupling problem. The case of the acoustic-acoustic coupling, which we are interested in, is 
modeled by the presence of spatially varying coefficients in the weak form of the equation 
m (see [2] for the linear and [3] and m for the nonlinear case) as follows: 

Find u such that 

f Q T ~ 2 k(x)u)ii(f> + Vu • V0 + b(x)( 1 - 6(x))Vu • V<£ 

+5(x)(5(x)|Vft| 9 ~ 1 Vh • V<p — ^yy-(f t ) 2< / ) } dxds = 0 

holds for all test functions <f> € X, 

with (it, h)|i=o = (it 0 ) and appropriately chosen test space X. In this model b denotes the 

quotient between the diffusivity and the bulk modulus, while the other coefficients retain 
their meaning. For brevity, we emphasized the space dependence of coefficients in (11.31) . 
while omitting space and time dependence of u and the test function in the notation. 

1.1. Some difficulties related to the model (11.31) . In |T7], the problem of optimizing the 
excitation part of the boundary in lithotripsy is considered, with the initial-boundary value 
problem for the Westervelt equation (11.11) as the optimization constraint. This situation 
arises when excitation and self-focusing of high intensity ultrasound is performed by a piezo- 
rnosaic. Compared to the the problem investigated there, the case of focusing by an acoustic 
lens presents us with several additional challenges. Not only a part of the domain boundary is 
optimized, but a subdomain which lies in the interior of the domain; this implies providing 
shape sensitivity analysis for an acoustic-acoustic interface problem. Insufficient spatial 
regularity of the primal (at most W l,q+l in space) and the adjoint state ( H 1 in space) on 
the whole domain does not allow for the shape derivative to be expressed in terms of the 
boundary integrals. However, it turns out that the state variable exhibits H 2 -regularity on 
each of the subdomains, provided that its gradient remains essentially bounded in space and 
time on the whole domain, which we will be able to use in our advantage. 

Working with the state equation also implies handling the nonlinear damping term of the 


(1.3) 
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(time-dependent) g-Laplace type, which is in itself a nontrivial task. Results on hyperbolic 
equations involving a g-Laplace damping term are sparse, and to the authors’ best knowledge 
do not include considerations of interface or shape optimization problems. 

1.2. Notational remark. To avoid confusion, it should be emphisized that we use a dot 
notation for time differentiation and t E [0, T] for the physical time variable, while r E M 
will be used for the artificial time variable to indicate varying domains. If 12+ is the initial 
shape of the lens, then 12+ iT will denote the perturbed lens obtained by moving points into 
the direction of some vector field h by some steplength r. 

1.3. Outline of the paper. The rest of the paper is organized as follows. We formulate 
the shape optimization problem in Section [2] and perform an analysis of the state and the 
adjoint problem. Section [3] deals with the question of existence of minimizers. In Section [4] 
we investigate the state equation on the domain with the perturbed lens. Some preliminary 
results needed to obtain the shape derivative of the cost functional are given in Section [5l 
Finally, Section [6] contains computation of the shape derivative. 


2. Shape optimization problem 

Let 12 C M d , d E {1, 2,3}, be a fixed bounded domain with strongly Lipschitz boundary <90, 

and 12+ a subdomain, representing the lens, such that 0_|_ C 11 and 12+ has strongly Lipschitz 

boundary <912+ = T. _ 

We denote by 12_ = 12 \ 12+ the part of the domain 

representing the fluid region. We then have <912_ = 

TU<912. 

n+, n_ will stand for the unit outer normals to lens 
12+ and fluid region 12_. Restrictions of a function 
v to 12 +i _ will be denoted by u+, V- and [a] := 
u+ — V- will denote the jump over T. 

Note that the assumptions on the regularity of the 
subdomains will eventually have to be strengthened 
to C 1,1 in order to express the shape derivative in 
terms of the boundary integrals over the interface 

r. 


CL 


Q 



Lens Q_|_ and fluid regions 


We consider the following optimization problem 


( 2 . 1 ) 


mm 

o + eo a d 

n€L 2 (Ox[0,T]) 


J(u, 12+) = 


mm 

fReCha 

n€L 2 (Ox[0,T]) 


(u — ad) 2 dx ds 


subject to the constraint (1.3), 


with the test space X = L 2 (0, T; 1L 7 0 1 ’ 9+1 (12)), and where the coefficients A, k, g, b and 5 are 
piecewise continuously differentiable and allowed to jump only over the interface <912+ = T, 


i.e. 


( 2 . 2 ) 


A, k, g,b,5 £ L°°( 12), 

Wi := w\n t € C' 1 (12j), for w € {b, g, A, <5, A;}, i € {+, -}. 


Ud € L 2 (0,T;L 2 (12)) is the desired acoustic pressure; O a d represents the set of admissible 
domains and is defined as follows: 


d a d = {12+ : 12+ C 12, 12+ is open and Lipschitz with uniform Lipschitz constant Lq}. 
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We assume that q > 1, q > d — 1 for the state and the adjoint problem to be well-posed; 
however we will need to strengthen this assumption later to q > 2 in order to prove certain 
properties needed for the characterization of the shape derivative. 

The strong form of the PDE constraint is given by 

r (1 - 2 k(x)u)u - div(^yVu) - div(6(x)((l - 5{x)) + 5(x)|Vu| 9_1 )Vu) 

_ 2k(x) ^2 j n |j 


(2.3) 


M x ) 

ui = 0 on T = 

1 du 


- «(*))$: + 6(x)J(x)|V»|«->^- 

u = 0 on dLl, 

^u,u)\ t =o = (u 0 ,ui). 


= 0 on T = <9fT 


For simplicity of exposition we consider homogeneous Dirichlet boundary conditions on the 
outer boundary <9fh We mention in passing that Neumann and absorbing boundary condi¬ 
tions as in mi, m could be easily incorporated here as well, with some additional terms 
involving the Neumann boundary excitation in the analysis of the state equation, but with 
no changes in the shape derivative itself, since the outer boundary is not subject to modifi¬ 
cations. 


2.1. Analysis of the state equation. Let us now add assumptions regarding the sign of 
the coefficients in the state problem: 

' X,k,g,b,S G L°°(ft), 

^ ( Wi := tclni G C' 1 (fl i ) for w € {b, g, A, 5, k}, i € {+, -}, 

w := \w\L°o(ty for w € {b, g, A, 6 , k}, 5 < 1, 

3g, b, 5 : A > A > 0, g > g > 0, b > b> 0, 5 > 5 > 0. 

We will utilize the following well-posedness result (cf. Theorem 2.3. and Corollary 4.1, [3]): 

Proposition 2.1. (Local well-posedness of the state problem) Let q>d — 1, q>l and the 
assumptions (|2.4I) hold. For any T > 0 there is a kt > 0 such that for all uq,u\ € W 0 1,9+1 (f2) 
with 

l u ili 2 (0) + I Vm oIl2(q) + |Vui|i 2(n) + |Vuol2 9 +i(n) + |Vui|^+i( n ) < k ? 

there exists a unique solution it G W C A = H 2 (0, T; L 2 (Ll)) n W 1,oo (0, T; W 0 1,,?+1 (f2)) of 

f| 1.3 h . and 

(2.5) W = {v € A : ||h|| i 2( 0i T;L 2 (a)) < ™ A ||Vu \\ l oo ( p , t - l ' 2 (Q)) < ™ 

^II^IU®+l(0,T;L®+l(n)) < M}, 

with 

^C^ r i, q +i LO o( KT + Ti+iM) < 1 , 

and fh sufficiently small. Furthermore, the following estimate holds: 

IMli°°(0,T;L°°(ft)) + INll 2 (0,T;L 2 (ST)) + ll^ 7 '“lli 2 (0,T;L 2 (n)) + II Vh| (o.T^+hH)) 

(2.6) + |H||oo( 0)T;i 2(Q)) + l|Vu||^oo( 0)T;i 2(Q)) + ||Vu||^<x.( 0)T;L 2(Q)) + ||V«||^to( 0iT;L 9+l(q)) 

< C(\ui\f2(ty + |Vito1 12 (q) + |Vui|| 2(f2) + |Viti| 

1®+ I (n))’ 
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where C depends on A, g, b, 6, k, and the norm x Li of the embedding Hq (D) ^a L 4 (fl). 
Since u — u,i € L 2 (0,T; L 2 (D)), the cost functional is well-defined. 

Degeneracy of the Westervelt equation is avoided here by employing the embedding 
W 0 1,9+1 (fi) ^A L°°(D) (since q + 1 > d), and the following estimate 

< C^, q+l L jVu(t)\ Lq+ i {Q) 

— ^VUL <?+1 roo + [ V'ii(s) ds |l/ 2 +l(Q) 
o j o 


- C w^ +1 , 


(|Vu 0 | L9+ i ( o) + (t q j" J^\Vu(y,s)\ q+1 dydsy +1 ^, 


which leads to the bound 

1 — oo < 1 — 2 ku < 1 + ao, 
( 2 -7) . . oM 


ao 2kC^i ,9+1 x,oo(l^ 7 “o|L«+ 1 (n) + T q+1 ll^ 7 '“|| L q+l(0 l T;I,9+ 1 (n)))- 


Due to the Sobolev embedding T / F 1 ’' ?+1 (D) ^A C 0,1 9+!(D), we know that u is Holder con¬ 
tinuous in space, i.e. u € C o,1 (0,T-,C°’ 1 9+1 (H)). 

2.2. Inequalities. Before proceeding further, let us recall several useful inequalities which 
will help us deal with the nonlinear damping term appearing in the state equation. From 
now on, we will always use C q to denote a generic positive constant depending only on q. 
At several instances, we will utilize the following representation formula for vectors x, y € 

(cf. Chapter 10, [23]): 

\ x \ q -l x -\y\ q -ly 

(2.8) r 1 r 1 

= (x-y) \y + v{x-y)\ q da + (q- 1) £(y + a(x - y),(x - y)da, 

Jo Jo 

where 

(2.9) £(x,y) = \x\ q ~ 3 (x-y)x 
as well as the inequality 

(2.10) | y + cr(x - y)\ q ~ l < \y\ q ~ l + |x| 9_1 , a € [0,1], 

For any y > 0, q > 0 and x,y € R d , the following inequality holds (cf. Lemma 5.1, [25]): 

(2.11) Wxr'x - \y\ q ~ l y\ <C q \x- y] 1 ^ (\y\ q ~ 1+ri + \x\ q ~ 1+q ). 

From here we also get 

(2.12) llxr 1 - I y\ q ~ l \ <C q \x- 2/| 1_r? (|2/| <? — 1+77 + \x\ q ~ 1+rt ), 
for 0 < y < 1, q > 1. 

We also recall that the following estimate holds for q > 1: 


(2.13) 


(M 9 1 x-\y\ q 1 y)-(x-y)> 2 1 q \x - 


y\ q+l > 0. 
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With the notation (12.91) . as a simple consequence of (12.111) . we have for vectors x, y, z, w and 
any 77 > 0 , q > 2 that 

I £(x,y) ~ C{z,w) | 

(2.14) = |(x • y){\x\ q ~ 3 x — \z\ q ~ 3 z) + \z\ q ~ 3 ((y — w) ■ x + w ■ (x — z))z\ 

< C q \x — z\ 1 ~ v {\x\ q ~ 3+ri + |z|' J,_3+?? )|x||y| + \z\ q ~ 2 (\y — w\\x\ + \w\\x — z\). 

We will also need Young’s inequality (see for instance Appendix B, m) in the form 

(2.15) \xy\ < e\x\ r + C{s,r)\y\ Vi' (e > 0, 1 < r < 00 ), 

r _ 1 

with C(e, r) = (r — l)r r ~ 1 e 1 ~ r . 


2.3. Analysis of the adjoint problem. Since the state equation contains a g-Laplace 
type damping term, its linearization provides a significant challenge. The main difficulty in 
considering the linearized q-Laplacian lies in the need for an essential bound on the gradient 
of the solution u] this is necessary for the linearized operator to be bounded. 

We will consider the adjoint problem with the assumption that the solution u of (11.31) 
exhibits the following regularity: 

(fti) «e tf 2 (0,r ; w 1 > oo (ft))n x, 


Note that the hypothesis is equivalent to assuming Lipschitz continuity of the acoustic pres¬ 
sure in space, i.e. u € H 2 ( 0,T;C 0 , 1 (ll))nI (cf. Chapter 2, Section 2.6.4, [5j)- Due 
to the embedding H 2 (0,T) W 1,oo (0,T ), if the hypothesis holds, we also know that 

u € W 1 ,oo (0, T; W 1 , 00 (f2)). 

The adjoint problem in weak form is then given as 


(2.16) 


7o r /«{*(! - 2 MK + iVp • VC - 6(1 - 6)Vp • VC 

< -b5(G u (S/p)y ■ \7(}dxds = f n j'(u)( dx ds 

holds for all test functions C € X', 


with (p,p)\t=T = (0,0), where j{u) = ( u — Ud ) 2 , X' = L 2 (0,T-, H 1 (Q)) and we have used the 
notation 


G U (Y) := |Vu| 9_1 Y + (q- 1)|Vu| 9 - 3 (Vti • Y)Xu 


for the linearized g-Laplace operator. The strong formulation of the adjoint problem then 
reads as: 


^(1 — 2 ku)p — div(^Vp) + div( 6 (l — 5)Vp) + div (6 6(G u (Vp))') 
= j'(u ) in 0 + U D_, 


(2.17) 


bl = ° on r, 

fi 9|L_6(1_ 5 )dp_ 

q on+ \ ' on.\- 


p = 0 on (912, 


bd ( G u {yp) ■ n + Y 


= 0 


on r, 


^(p,p)\t=T = (0,0). 


To obtain a forward problem, let t 6 [0, T] and p(t) := p(T — t) (cf. Lemma 3.17, [30]): then 
(p,p)\t=o = ( 0 , 0 ), and we have the following problem for p: 
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(2.18) 


7o T /n{i(l - 2 ku)k + JVp • VC + 6(1 - 6)VP • VC 
+6<5(|V , ii| 9_1 Vp - (| Vul^yVp) ■ VC 
-(q - l)66((|V7^- 3 )-(Vp • Vu)V« - \Xu\ q ~ 3 (Xp ■ Vu)Vu 
+ |V«| <? ~ 3 (Vp • Vii)Vu + |Vu| 9_3 (Vp ■ Vuj'Vil ^ • Vc| cfcrds 

= fo T J n j'(u)(dxds, 

holds for all test functions C £ X'. 


Proposition 2.2. (Local well-posedness of the adjoint problem) Let q > 1, q > d — 1, 
assumptions (El on coefficients and hypothesis fUi) hold. For sufficiently small m, final 
time T > 0, and llViiH^^g T-L°°<n)) I|V'6 ||l oo (o,t ; l oo (Q))i there exists a unique weak solution 
p£ X = C 1 (0, T-L 2 {n))nH 1 (0,T; H 3 {n)) of d2Tffl . 


Proof. The well-posedness of the adjoint problem can be obtained through standard Galerk- 
ing approximation in space, energy estimates and weak limits (cf. Section 7.2, [JJ). We will 
focus here on obtaining the crucial energy estimate. Testing (12.181) with ((a) = p(cr)x[o,t) 
(first in a discretized setting and then via weak limit transfered to the continuous one) yields 


(2.19) 

where 


1 

2 L 

< - 


/ (1 — 2ku(a))(p(a)Y dx + 

In Jn 6 


11 Vp(cr)| 2 dx 


[ [ 6(1 -5)\Xp\ 2 dxds 
Jo Jn 

n \-u{p) 2 dx ds + f f j'(u)pdxds+ f f bdp\Xp\\Xp\dxds, 

! * Jo Jn Jo Jn 


ip := KlVtil 9 - 1 )'! + (q- l)|(|V?i| 9_3 )'||Vu| 2 + 2 (q - l)|Vh| 9 - 2 |Vii| < C g |Vii| 9 - 2 | Vii|, 


and we have used that Jo Jo ^ l^l 9 *| Vp| 2 dx ds > 0 and Jq fn(q — l)65|Vu| 9 3 (Vp • 
Vh)(Vri • Vj5) > 0. By taking essential supremum with respect to f £ [0,T] in (12.191) 
and employing 

[ [ M|Vp||Vp| dx ds 

Jo Jn 

— C g ||Vu||^ oo ^ 0>T . LOO( . n ^||Vii|| L 2( 0jT . I/ oo(Q)) (^7ll^ 7 p|||oo(o,T ; i, 2 (a)) + £ II^ 7 pIIl 2 (o,t ; l 2 (Q))) > 

we get the energy estimate for p 


/I k \ 

( 4 ^ _ fl °) ~ ^’^'( ( ^H 1 ,L 4 ) 2 |I^IIl oo ( 0 ,T;L 2 (Q))J IIpIIl°°( 0 ,T;L 2 (Q)) 

+ (i^ ~~ ^^C , <?ll^ 7 ^llL° o (0,T;L°°(n))ll^ 7 ^Hi 2 (0,T;L oo (n))) II^PllI °°(0,T;L 2 (O)) 

+ ( 2 -^ “ 7 ~~ ^(C , /|i,L 4 ) 2 |I^IIl oo ( 0 ,T;L 2 (Q)) 

— £b5C q \\Vu\\ g Lx ^ 0T . Lx ^\\Vu\\ L 2^ T . La o^^j ||Vp|||2( 0i T;L 2 (n)) 
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< 


/ [ j'(u)pdxds < 2\\u-u d \\ L 2 {0tT . L 2 m \\p\\ L 2 {0tT . L 2 m , 

Jo Jn 


with 1 — ao > 0, defined as in (12.71) . bounding away from zero factor 1 — 2 ku. This estimate 
holds under additional assumptions on smallness of T, fh and || Vu |||^ 0 T-L°°(n)) ll^“IU 2 ( 0 ,T;L oo (n)) 


( 2 . 20 ) 


Tj(C^^) 2 Cprn<^(l-a 0 ), 

-b5C q ||V«||« oo(0>r;2 . oo(n)) \\Vu\\ L > { 0 ,T;L~ m < =■ 


1 9-2 

uu^y a || v u,| 

£ 

\{C%^) 2 Cpfh<\b(l-S), 


' 2 J 

and some sufficiently small e > 0. Here we have employed Poincare’s inequality for functions 
in Hq(£1) (see Theorem 2.17, [23] ): 

INlL°°(o,T;Z, 2 (n)) < Cp||Vu|| L oc( 0 , T;L 2 (n)) < Cprh , u E W, 

with Cp = Cp(H) > 0. Moreover, we have 

(p(*)» 0 i, 2 (n) 


IpWIh-ho) — 


sup 

C6Hi(n)\{o} 

-h 


( 2k 1 

i -r U P( ~ -Vp ■ VC + 6(1 - <5)Vp • VC 
1 A Q 


- 6<5(|Vu| 9_1 Vp - dVul 9 " 1 )'Vp) • VC 




+ (9- 


9—3y 


(Vp • Vu)Vu — |V-u | 9 3 (Vp-Vit)Vu 


+ |Vti| 9_3 (Vp • Vii)Vu + |Vu|' ?_3 (Vp • Vu)Vii) • VC + /(u)C} d®. 


This further implies 

|p(<)lff-i(n) < c '(l u ( i )U°°(n)|p(0l^- 1 (n) + |Vp(t)| L 2 (n) |Vu(t)|^“ 2 (n) |Vu(t)| L oo( n) 

+ IVpWIls^IVu^)!^^ + I u(t) - u d {t) | L 2 (n) ), 
with C = C(k, A, g, 6 , 5, 5 ) > 0. By squaring and integrating over (0, T), we can achieve that 

p : €i 2 (o,r ; r 1 (n)). 

According to Theorem 3, Section 5.9, [7], since p E T 2 (0, T; Hq(Q)) andp E L 2 (0, T ; i7 _1 (fl)), 
it follows that p E (7(0, T; L 2 (H)). The statement then comes from reversing the time trans¬ 
formation. □ 


3. Existence of optimal shapes 


From now on, for simplicity of exposition, we assume that all the coefficients in the state 
equation are piecewise constants, i.e. 


(3.1) 


' A ,k,g,b,5e L°°(n), 

< Wi := rein; is constant, for w E ( 6 , g, A, 5, k}, i E {+, —}, 
Wi > 0 for w E {b, g, A}, Si E (0,1), ki E M. 


Note that now u = min{|cj + |, |w_|}, u = max{|o; + |, |w_|}, where oj E {b, g, A, <5, k}. 

We turn next to the question of existence of minimizers for the shape optimization problem 
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m, with the coefficients satisfying assumptions (13.111 . The main idea of the proof is to 
employ the Bolzano-Weierstrass theorem for continuous functions on compact sets (for this 
approach see, for instance, Section 3, [9] and Section 3.2, m- 
We recall the following compactness result (cf. Theorem 2, [9j): 

Theorem 3.1. Let 12+ be a sequence in O a d. Then there exists 12+ £ O a d and a subsequence 
12+ fc which converges to 12* in the sense of Hausdorff, and in the sense of characteristic 
functions. Additionally, 12)J) fc and <912+ fc converge in the sense of Hausdorff towards 12^_ and 
<912+, respectively. 

Here the set of characteristic functions is defined as 

CharfTl) = {xp+ : 12+ C 12 is measurable A xn+(l — xn+) = 0 a.e. in 12} 

and convergence on this set first of all means pointwise almost everywhere convergence of 
functions, but due to Lebesgue’s Dominated Convergence Theorem also convergence in L p (12) 
for any p € [ 1 , oo). 

We first establish continuity of u(12+) with respect to the shape of the lens 12+. 

Proposition 3.1. Let q > 1, q > d— 1 and the assumptions (ED on the coefficients in (ED 
hold. Then the mapping xq+ 11 ( 12 +) is continuous from the set of characteristic functions 
Charifl) to W ^(O, T; L 2 (12)) n W^ +1 ( 0, T; W 0 1 ,,+ 1 (12)). 

If additionally there exists e > 0 such that the solution u ** of fll.3|> with 12+ = 12+ satisfies 

(3.2) lMlllU 1 .5+ 1 (0,T;IU 1 '9+ 1 + e (T2)) < C) 

where C depends only on 12 and the final time T, then the mapping xu+ | -A ii(12+) is even 
Holder continuous with exponent | at x^ + from the set of characteristic functions Char(Ll ) 
mL r2 ( 12), r 2 = (1+g ^ 1+£) to W 1 ,oo (0, T; L 2 (12)) n W 1 , 9 + 1 (0, T; W 0 1 , 9 + 1 (12)). 

Li 

Proof. Let 12+ be a sequence converging to 12+ in the sense of characteristic functions. By 
subtracting the weak forms for u n and vfi , corresponding to the domains with the lens regions 
12 + and 12 + respectively, we get 

pT p I o k- 

V / / {t-(1 - 2kiU n )(u n - - ^-{u 71 - - V(u n - u^) ■ Vfii 

Jn Xi Xl Qi 

+ bi( 1 - 6i)V(u n - it # ) • Vfii + Mi(|Vu n | 9_1 Vu n - |V« Jt |« _1 V* # ) • V0 
2k- 

- - {u\f)4>}x n? dxds 
Ai 

rT r \ 1 

= / / { — (1 — 2 kiU^)i$(j)-\ - \7u^ • S/fi + bfil — <5j)Vu^ ■ V(j> 

. c r, i do Jn *i Qi 

2k- 

+ bi5i\Vu‘^\ q ~ 1 X7v} ■ Vfi - (xn? - x n »)dxds, 

Ai t 

with 12” = 12 \ 12+, 12L = 12 \ 12+. Testing with <ft = ii n — and employing inequality (12.131) 
for the difference of the g-Laplace terms then yields 

— (1 — 2k\\u ||x,<x> (o,T; ;L°° (f2))) 11^ — ^ llL°°(0,T;L 2 (n)) + 7 ^\\^( u n ~ u m) \\l 2 (0,T;L 2 (Q.)) 
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+ I>(1 — ^)||V(li rl — Um) ||L 2 (0,T;L 2 (f2)) + ^ 9 MI|V(u — tl^) IIl<j+1(o,T;L9+ 1 (Q)) 
2k 

< ~dl^ tl ||L 2 (0,T;L°°(n))ll^ n - ^ # lli°°(0,T;L 2 (n)) 

+ (C// 0 1 ,L4) 2 ||n rl - ^11 L°°(0 ,T;Hq (Q)) II^ IIL 2 (0,T;L 2 (fi)) II~ ^\\L 2 (0,T]H^(CI)) 

+ V I [ [ {-^-(1 — 2kiU^)u\u n — vJ) + — SIvJ ■ X7(u n — viJ) 

.Jr 1 I Jo ./n l A Bi 

+ bi( 1 - <5i)Vii # • V(u n - u # ) + 6 i <S i |V« ll | 9_1 Vu lt • V(ti n - u # ) 

- ^(J) 2 (u n -u^(xn^ ~ X^)dxds . 


By employing the Sobolev embedding W 0 1,9+1 (fi) L°°(Q): 

II^IIl 2 (0,T;L°°(T2)) < < - ! Wo’ q+1 ,L°°^^L 2 {0,T-,Wg’ q+1 (n)) 

and then utilizing (12.61) to estimate the terms 11 ^ 11 ^, 2 ( 0 T-w 1 ’ 9+1 (f2)) anc ^ 11^1 lL 2 (o,T;L 2 (n))) we 
conclude that for sufficiently small initial data the two first terms on the right hand side can 
be absorbed by the appropriate terms on the left hand side. For the remaining terms, we 
employ the following estimates: 

/ / {tU ~ 2k i 'J)u\u n - u # )(xQ" - X a t)dxds 

Jo J!i u 

< (1 + 2 fc||u t) ||L°°( 0 ,T;L°°(n)))ll^ ti ||L 2 ( 0 ,r;L 2 (n))ll 'U'n &m IU 2 (0,T;L9+ 1 (n))IXfi? — X n #Un(0)> 

together with 




n 


« # xxn? 


Xq| ) dx ds 


2k 

< — 

“ A 


\W 


|| 2 

llL tx> (0,T;Z,°°(f2)) 


~ ^ ll-L 2 ( 0 ,T;L 2 (f 2 )) |Xfl™ 


X n #lL 2 (f2), 


and the estimate 
T 

I [ -VJ ■ V(u n - u^)ixn n - Xa) dxds 

Jo Jn 0 1 '■ 

< -|| Vw# llL 2 (0,T;L 2 (n))l|V(u n - U # )||L 2 (0,T;L'?+l(O))|Xn7 ~ XQ#lL r (f2) 

Q i 

< -To ||Vu # || L 2 (0iT;i 2( n) )||V(u n — 'h tt )||L9+i(o,T;L'?+i(Q))|xn^ - Xn#U ri (fi)’ 
0 * 


with n = 2 ( ff } . An analogous estimate to the last one can be derived for the 6j(l — <5j)-term 
where Vu^ is replaced by ViiA Thus all these terms on the right hand side tend to zero as 
n —>• oo. Finally, by Lebesgue’s Dominated Convergence Theorem also the integral 

rT r 

biSilvJr'Vu* ■ V(ii n - «#)(xn- - Xn t)dxds, 


j f 

Jo Jn 


whose integrand due to the factor (xnv- — Xq#) tends to zero pointwise a.e. and is bounded by 

the integrable function WUVudiyi^.^+i+e^))(|| Vti n || L9+ i( 0!T;Z/9 +i ( Q )) + ||Vu tt ||L9+i(o,T;L-J+i(f2)))xn, 
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goes to zero as n -A oo, which proves the assertion. 

If we additionally assume that the solution of the state problem exhibits a slightly higher 
regularity in space than IT 1 , 9 + 1 (D), namely that u € VL 1 , ' ?+ 1 (0, T; W 1 , 9 + 1 +e (D)) for some 
e > 0 , we can enhance the estimate on the latter term to 

[ T [ -V(u n -ul){xm-X n t)dxds 

Jo Jn 1 

— bSWVu^ ll^g+l(o ! T;£ / 9+ 1 +e(f2)) ll^(^ — lli , 3 , +l(0,T;L | J+ 1 (Q)) IXfi™ — Xq# \L r 2(p), 

with r 2 = ^ 1+<? )( 1+g ) . Altogether, we can then conclude that 

Q£ 

(1 — 2k\\u ||l / oo (0,T;L oo (f2))) IIh — h^|| i oo(0 i T;I, 2 (n)) + ll^( u — n **)llL 2 (0,T;L 2 (n)) 

+ ll^(h n — U^) Hl 2( 0! T ; L2(Q)) + l|V(h n — ^)||^+l( 0 ,T;L9+l(n)) 

<C((l + 2k\\u^ ||Lo°( 0i T ; L°°(n)))ll'“ tt H l 2 (o,t ; L 2 (n))||h n - ^||L 2 ( 0 ,T ; T?+ 1 (fi))lxn^ — x^duifri) 

+ ll^ 7 ' utt ||L 2 (0,T;L 2 (a))l|V(u n - U tt )|| M+ 1 (o,T;L9+i(n))|xn^ — Xq# |z,n(n) 

+ ll^ 7 '“^llL 2 (0,T;L 2 (a)) l|V(^ n ~ IIl , J+ 1 (0,T;L5+ 1 (O)) |Xf2" “ Xq# |z,n(f2) 

+ || ||^ 9+ ij.o i 2 l ; £ / 9+i+e(Q)) l|V(lt — ||L9+l(0,T;L9+i(f2)) IXfi” — XQ#jL r 2(n) 

+ ll u ”ll|^(0,T;L°°(f2))ll 7 hi “ ^m||L 2 (0,T ; L 2 (n))IXn^ — Xn# Il 2 (Q)), 

for sufficiently large C > 0 independent of n, which implies Holder continuity of the mapping 

XH+ u(D + ). □ 

Now let H” be a minimizing sequence for (|2.1D . Due to Theorem l3.ll there exists a subse¬ 
quence, which for brevity we still denote D”, that converges to some Di. € O. By extracting 
another sequence, we may as well assume that Dll -A D* in the sense of characteristic func¬ 
tions. 

Let us denote by u n the weak solution corresponding to the domain where the lens region 
is given by D”. We know then that u n satisfies the estimate (12.61) . where u is interchanged 
with u n . This means that we may extract a subsequence, again denoted u n , such that 

u n u * in H 2 (0, T ; L 2 (D)), 

Vu n — 1 Vu* in L 9 + 1 (0,T;L 9 + 1 (D)). 

Due to the embedding H 2 (0,T) ^A C 1 (0,T), this further implies that u(t) —*■ u*(t) and 
u(t) —- u*(t) in L 2 (D) for all t € [0, T], Therefore, (u*,u*)|t=o = (uo,ui). It remains to 
show that u* solves the state problem on the domain whose lens region is given by Dlj_. 

Proposition 3.2. Let the assumptions of Proposition \3.1\ be satisfied. Let D” be a min¬ 
imizing sequence for the shape optimization problem m and let D* be an accumulation 
point of this sequence in accordance with Theorem 1,9.11 Then the sequence u n correspond¬ 
ing to the domain with the lens region D” converges strongly to u * in VL 1 ’ oo (0, T; L 2 (D)) n 
W 1 , 9 + 1 (0, T; Wq 1 , 9 + 1 (D)) ; where u * is the solution of (11.31) in the domain whose lens region 
is given by D+. 

Proof. In order to see that u * is the weak solution of the state problem in the domain where 
the lens region is given by Qlj_, note that due to the fact that u n solves (11.31) with lens region 
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0” and using integration by parts in the first term on the right hand side 

V f f { T-(! - 2 k iU *)u*(j) + —Vu* • V0 + bAl - 5 t )Vu* • V0 
v-7^ i Jo Jn h Qi 

2k- 

+ bidi\Vu*\ q ~ 1 \/u * • V0--2-(h*) 2 ((>}xn* dx ds 

^i % 

rT r 1 r)U, 1 

. c 77 - ' , jo in Aj A* Qi 

(3.3) + 6i(l - 5*)V(h* - u n ) • V0 + 6^i(|Vu*| 9_1 Vu* - iVhT^Vfd 1 ) • Vc j> 

2k- 

- fTl(u* - h n )(h* + u n )cf>} xn* dxds 

1 

+ V f [ {-^-(1 — 2kiU n )u n 4> + —Vtt n • V0 + 6j(l — (5j)V'ii ri • V(/> 

^V° Jn A, ft 

2k- 

+ b i 5 i \Vu n \ q ~ i Vu n -V(t)- —2-(?t n ) 2 <?!>}(xrR ~Xw)dxds, 

Aj 

for all 4> € C°°(0, T; Cg 0 ^)), <)>(T) = 0. The difference of the g-Laplace terms on the right 
hand side can be estimated with the help of inequality (|2. IIP (with r] = 0) as follows: 

[ [ 6 i <y i (|Vu*| 9_1 V«* - |Vu”| 9 " 1 Vh n ) • Vfxn* dxds 
io Jn 

<C q bd f I \V{u*+ \Vu n \ q - 1 )\V0\\xa+\dxds 
J o Jn 

<Cqb5\\\7(il —u ) ||l5'+1(0,T;L9+ 1 (Q)) II L‘2+ 1 (0,T;L | 3+ 1 (Q)) 

+ ll V?i ' l |ll,<J+l(0,r;L'?+l(n))) ll V< ^lll,9+l(0,T;L9+l(n))- 

The remaining terms can be estimated analogously to the estimates in the proof of Propo¬ 
sition EU from which it then follows that the right hand side in (13.31) tends to zero as 
n —> oo. □ 

Theorem 3.2. The shape optimization problem (E3D has a solution. 

Proof. Let us define the reduced cost functional J : O a( j —>• M: 

J(n+) = J(u(n+),n + ). 

Let as n -2- oo. It can be shown (cf. Lemma 3.3, [2Tj) that 

| J(u*, 12* ) - J(u n , fi™)| < ||?X* — l/ ?! '||L 2 ( 0 ,r;L 2 (n))ll' lt * + — 2Ud|| i 2( 0j T;L 2 (fi))- 

Since the term \\u* + u n — 2ud\\L 2 (o,T;L 2 (n)) is uniformly bounded due to (12.61) . by employing 
Proposition 13.11 we achieve that the right hand side tends to zero as n -2- oo. Therefore the 
cost functional J is continuous on O a d- According to the Bolzano-Weierstrass theorem, since 
O a d is compact, J attains a global minimum on O a d- □ 
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4. State equation on the domain with perturbed lens 

4.1. Method of mappings. The approach to computing the shape derivative that we will 
take follows the general framework given in p3j and its extension to a time dependent setting 
given in m- One of its main ingredients is the mapping method, originally introduced by 
Murat and Simon in [26] . which we briefly recall. 

We introduce a fixed vector held h G C' 1,1 (fi, R d ) with h \qq = 0, and a family of transfor¬ 
mations F t : Q i—)• M. d 

F r = id + rh, for t G M. 

There exists to > 0, such that for |r| < to, F t is a C 1,1 -diffeomorphism (cf. [5] ) . If the 
perturbed lens is given by fl +iT = F r (Q + ) and T r = F r (T), then it follows that T r is 
strongly Lipschitz continuous (cf. Theorem 4.1, [12]). 

The Eulerian derivative of J at f2+ in the direction of the vector held h is defined as 

dJ(u , f 1+)h = lim —(J(u T , f2+ r ) — J(u, S4+)), 

r—>0 T 

where u T satishes the state equation on the perturbed domain fl r . The functional J is said 
to be shape differentiable at fl_|_ if dJ(u,Ll + )h exists for all h G C' 1,1 (fi,M d ) and defines a 
continuous linear functional on C' 1,1 (f2,M d ). 

We introduce the following notation: 

(4.1) I T = det {DF T ), A r = ( DF r )~ T . 

where DF r is the Jacobian of the transformation F r . 


Lemma 4.1. [14] For sufficiently small tq > 0, mapping F r has the following properties: 


o irbdxd 


t i— >• F t g C 1 (—To,T 0 ;C' 1 (n,M d )) t i-a A r G C(—r 0 , to; (7(11, 

T eA F t g C(— tq, tq; C 1,1 ^, M d )) tg/ t G t 0 ; C(ty) 


)) 


t i-a FF 1 G C(— tq, tq; C 1 (n,M d )) F 0 = id 


d ^ , 

— E r=0 = h 

dr 


d_ 

dr 


DF r | t= o = Dh 


±F 

dr 

d 


-ii 


=n = —h 


I T— 0 


-^-I t \t=o = di vh. 
dr 


ir DF "\ - = 5F( A 

f-A T | T=0 = -(Dhf. 

dr 


) T |r=o = -Dh 


As a consequence of Lemma 14.11 there exist «o,Q:i >0 such that 

(4.2) 0 < ao < I T (x) < a i, for x G 12, t G [—To, to]. 

Furthermore, there exist /3i, /?2 >0 such that 

(4.3) |-4t|l°°(q) < 01, l^r 1 |l°°(Q) < 02, for t G [-t 0 ,t 0 ]. 

We will often employ the following lemma which gives us the formula for the transformation 
of domain integrals: 

Lemma 4.2. [29] Let <p T G L 1 (U r ), then ip T o F T G L 1 (fl) and 
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Let us introduce, for some Lipschitz domain € O a d, the operator E(-,Q + ) : W —>• X * 

by 

(E(u,n+),<f>)x ir % = V / / {t-(1 - 2kiu)u<j) + —Vu • V</> + 6j(l - dj)Vu • V0 

’ iG{+ _ } do dn. ( A* ft 

2k- 

+ 6j(5j|V-u| 9 _ 1 V'u • Vd>-- 2 -(u) 2 d>} dxds. 

A^ 

Fix r € (—To, To). Proposition 12.11 guarantees that E(u,Q +tT ) = 0 has a unique solution, 
which we denote u T : D —>• R. We transport u r back to the domain with the fixed lens f2+ 
by defining u T : fl —> R as 

u T = u T o F t . 


Differentiability of h implies that u T € X. We further have 
(4.4) 

{E(u t , n+, T ), db -} X*,X t 

rT 


y / / {—(i - 2/cju T )u T ^ T h— vitr-v^r 
1/0 1 A * & 

2/c ■ 'i 

+ 6*(1 — 5i)Vu T ■ X4 > t + 6j<5j|Vu T | 9 ^ 1 Vu T • V<^ r - —^-(u T ) 2 4> T > dx r ds 

A; J 

V f / f — (1 _ 2fciU T )ii r (/. T +-A r Vu T • A T V^ T + 6 i (l-<5 i )^ T Vu T • A r V0 T 
, , do dn^ 1 A* ft 




2k T 'i 

+ Mi|A r Vu T | 9 - 1 A T Vu T • A r Vd) T - — (u r )V [ I T dx ds 
= {E(u t ,t),4> t )x^x, 

for any cj) T € X T = X, with D_ iT = Q \ f2 +iT , where we have used Lemma 14.21 and the 
fact that Vu T = A r Xu T o i^T 1 . Therefore, for sufficiently small |r|, u T uniquely satisfies an 
equation on the domain with the fixed lens: 

(4.5) E(u T ,t) = 0. 

Since Fq = id, we have that u (> = u and E(u, 0) = E(u, Q + ). 


4.2. Continuity of the state with respect to domain perturbations. We will now 
focus our attention on the speed of convergence of u T to u as r —» 0 and prove two properties 
which together correspond to hypothesis (H 2 ) in [T4j (see also Proposition 3.1 in m- 
We begin with the question of uniform boundedness of u r with respect to r 6 (—To, To). 
Since 

11^ llz,°°(0,T;L°°(r2)) = I \ u t ° Fr ||L°°(0,T;L°°(Q)) 

< 11 11 L°° ( 0 ,T;L°° (H))! 

we also have that 

(4.6) 1 — ao < || 1 — 2ku T \\ L oo( 0tT . L oo(ty' ) < 1 + ao, 

where ao is defined as in (12.71) . 
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Proposition 4.1. Let q>l, q > d — 1 and assumptions EB hold. Solutions u T of (14.51) 
are uniformly bounded in 1P 1 , 00 (0, T; L 2 (Ll)) n W 1,q+l (f), T; VFg 1 , 9 + 1 (0)) for t e (—To, to). 

Proof. It can be shown (cf. Lemma 3.3, [13]) that 

1 + To|-D/l|£oo(f2) 


Vu T 


\L<i+ 1 (0,T;Li+ 1 (Q)) < 


a, 


l/(9+l) 

0 


Vu T 


li‘J+ 1 (0,T;L9+ 1 (n)). 


with «o as hi (14.21) This implies that, for u T £W, we can estimate 


(4.7) 


I Vti T 


Il 9 + 1 (0,T;L'?+ 1 (Q)) < 


1 + To|L>/i|loo(q) 


a. 


l/(9+l) 

0 


M. 


Testing (I4.5P with E as in (14.41) with <j) T = u T X[o,t) € L 2 (0, T; 14 / 0 1 ’ <?+1 (H)) yields 
1 


1 o i* i r r 1 o i* 

- -• T\t„-.T \2 T j_ + -\A T Vu T \ 2 I T dx 

j0 2 L J n Q Jo 


A 


(1 — 2 ku T ){u T YI T dx 


+ f [ 6(1 — <5)|A r Vu r | 2 I T dx ds + f f 65|H T Vu r |' ?+1 / T <ix(is 

Jo Jn Jo Jn 

= j f — (u T ) 3 I r dxds. 

Jo Jn A 

By taking the supremum over t € [0, T] and by utilizing the uniform boundedness properties 
621), ESP and estimates (|4.6j) and (14.Tp , we find that 

1 «0 /-. \ 11 • r 112 111 ||_ r ,|2 

4X _ ao) " U IIl°°( 0, T;L 2 (ft)) 4" VU " L °° (0,T;L 2 (ft)) 


+ 2“^ ^ ^2«0||V^||i 2(0iT;i 2 (n)) + 2^^+I a °ll V “ T |lLtf 1 (0,T;L?+ 1 (n)) 

~k 1 

— Ol ||'6 T ||L 2 (0,T;L°°(Q)) V^T||tl' |lL«=(0,T;L 2 (n)) + ^(1 + a o)«l Wl 1L 2 (p) 


< 


(^■8) + 9 ^ a i|Vuoli2( n ) 

- ^ ai ^ 1 .«+ i jioo r9+1 -i/(<Ji) ( ) ^ll' uT lli oo (0,T;L 2 (n)) + + a o)ai|ni|| 2(a ) 

+ —ai|Vu 0 || 2 ( n ). 

From here, for sufficiently small M and T, we can achieve that the first term on the right 
hand side gets absorbed by the appropriate term on the left hand side. Thus we have uniform 
boundedness of u r in W^°°{0, T; L 2 (Q)) n W , 1 >« + 1 (0,T; W 0 1 ,,+ 1 (n)), |t| < t 0 . □ 

The Holder continuity of u with respect to domain perturbations is established by our next 
theorem. 

Proposition 4.2. Let q > 1, q> d— 1 and let assumptions EH) hold. Then 

lim-(lK — ^ll|°°(0,T;L 2 (O)) + _ ' U ')IIl°°(0,T;L 2 (O)) 


(4.9) 


+ ||V(u T f0llL2(o )T; £ 2 (Q)) + ||V(u T u) ||^t+l(o,T;L 8 + 1 (f 2 ))) ~ ® 
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Proof. Note that the difference v T = u — u r satisfies 

rT 


y 


/ o Jn, 


1 2k 1 

" - 2 kiUi)vJfair ~ —^-vfUifilr H- Vvf ■ V(/>i 

■i Ai Qi 




(4.10) 


* 6 {+ 

+ bi( 1 - 8i)VvJ • V& + - iVufl^Vuf) • V& 

2 h 


\ \ 


dx ds 


(iii + uf)vj <fo} 

= (/+></’+)X*,X+ + (/-) 0-)x*,X_> 
for all 4> € X, with the two terms on the right hand side given by 
(/*) Wx*,Xi 

= / / |p(4 - 1)(1 - 2kiUi)ui<l>i + — ((/ r - l)H T Vu[ • 

+ (H T - /)V< • ArVfa + V< • (Ar - /)V0i) 

+ 6 i(l - <50 ((/ r - l)A T Vii[ • ArV<k + (A t - J)Vu[ • H r V<^ + Vu[ • (A T - /)V&) 

+ Mi((|A r V<| 9_1 A r Vu[ - |V«[| ,_1 V«J‘) • V& + |A r Vul’| 9_1 A r Vu[ • (A r - J)V& 

+ (I T - 1 )| ArWuU^ArWuJ ■ A T V(/>iJ ~ (I T ~ l)^(u[) 2 ^} dxds, 

where Xi = L 2 (0, T; W 1,q+l (fAi)). Since 4>i = ul"X[o,t) £ -X’, we can use it as test functions in 
(J4.10I1 . which together with the uniform boundedness of I T results in 

{ \ a ° [ T-(! — 2kiuJ)(vl) 2 dx + ^ j — \VvJ\ 2 dx] 

_ r | l */ ilj i </ Qi ft 


*e{+,-} 

rt 


JO 


+ [ f b i {l-5 i )\Xv T i \ 2 dxds + [ [ 60i2 1 ~' ? |Vu[| 9+1 dxds) 

Jo Jn t Jo Jn, -* 


0 -/D; 


* E { 

ie{+ 

■yfli 


■T II 2 


«l|Fi llL 2 (0,T;L°°(n i ))ll ?; i Hl^O.T;!, 2 ^)) 

+ 2(C“\ ^ 4 ) 2 Cll ||ifi ||x ,2 ( 0 ,T ; L 2 (r 2 i)) ll^r III / 00 (0,T-yH 1 (fi;)) ll^i IIl 2 ( 0 ,Tj //" 1 (fi;)) 


+ 2 (ll'“IllL 2 (0,T;L°°(n i )) + 11 11 A 2 (0,T; A°° (fi;)) ) 11 11 A 00 (0,T; A 2 (SI;)) 

Here C^\ stands for the norm of the embedding i5 1 (Hj) L 4 (flj) and we have also uti¬ 
lized the fact that (|Viij|« _1 Vuj - \X uf) ■ X{ui - uj) > 2 1 ~ q |V(?i; - itj)\ q+1 , which 

follows from (12.1311 . 




We can employ 

(4-11) ll^rilA 2 (0,T;L°°(r2i)) S J ' 2 W+D G '"‘ li9+1 Loa ||lt- || M + 1 ( 0 ,T; VE^+l^)), 

and the same inequality with u t instead of uj, as well as 


< x 2(9+1) c?. z * 
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(since v[\t=o = ^IIt=o = 0 ) to conclude that, because ||h)"||£ / 2 +i(o,T;W 1 ’'?+ 1 (n i )) un if° rm ly 
bounded for r € (—to,to), for sufficiently small m and final time T it holds 


■•T II 2 


(4.12) 


(II t ,T II L°° (0,T ;L 2 (0)) + ll^'Hli °°(0,T;L 2 (n)) + ||Vh T , 

lL 2 (0,T;L 2 (a)) 

*£{+,-} 

+ l|Vhl 

'L+ 1 (0,T;IA+ 1 (n))) 

< C(\(f+, v\)x* x + \ + \(f-,V-)x*,X_ l)» 


for some sufficiently large C > 0 which does not depend on t. By employing the inequality 
( 12 . 111 ) (with 7] = 0 ) we obtain 


(4.13) 


KlArVuJl^ArVuJ - |Vu[| 9_1 Vn[)Vhj| 

< c q \{A t - I)VuJ I (iVhir - 1 + I (A, - JjVtkTr 1 ) |Vhj 

We can therefore estimate the terms on the right hand side in (14.121) as 


(4.14) 


l(/i> t’i ) X* x i I 

< ^l^r - l|L°°(Qi)(l + ao)||ili||L 2 (0,T;L 2 (n i ))l|t ; i r ||L 2 (0,T;L 2 (O i )) 

+ (I A- ~ l|L°°(nj)^i + \A r — /|l«>(q.)( 1 + ^i))(~l|V«i"||L 2 (o,T ; L2(n i )) 

+ 6(1 - S)\\VuJ\\ L 2( 0 jT . L 2(q.^ ||Vh[|| L 2 ( 0 ,T ; L 2 (n i )) 

+ &<5^| 4l r — /|L«>(n i )(l + | A t — ^"|^oo(q.)) + \A t — 

+ l-^r — l|L°°(a i )/5i )l|Vhj ||| 9 +l( 0) T;L9+ 1 (n i ))ll^i ||L9+ 1 (0,T;L?+ 1 (f2 i )) 


ir il 2k 

+ Mr — 1 TOO in.-IT" 




(Hi) — W_fl+L 4 ' 11+ llz, 2 (0,T;.ff 1 (n i ))ll + llL°°(0,T;L 2 (£+))’ 


where a o is defined as in (|2.7I) . By inserting this into the estimate (14.121) and by employing 
the uniform boundedness result from Proposition su properties of the mapping F r from 
Lemma [4~D and Young’s inequality, we can conclude that 


^“5) ^ (ll^ T llz/°°(0,T;L 2 (f2)) + ll^ 7t ’ r |li o °(0,T;L 2 (n)) + Hl 2 (0,T;L 2 (T2)) 

(4.15) * 6 {+,-} 

+ ll^'W IIli+ 1 (0,T;L9+1(Q))) = 0‘ 

In oder words, we know that lim u T = u in X. To obtain the statement of the Proposition, 

T —>-0 

we divide (J4.12I) by t, and then it remains to show that 

I (/+> ^+) X*,X+ I + \(f-’V-)x* ,xA 

lim-±—-—-= 0 . 

t->0 T 

This now follows from the estimate (14.141) . Lemma 14.11 Proposition 14.11 and (14.151) . □ 


If we assume higher spatial regularity of u, we can even obtain Lipschitz continuity with 
respect to domain perturbations. 
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Proposition 4.3. Let q > 1, q > d — 1 and let assumptions hold. Assume that the 

solution u of satisfies 

(4-16) l|Vit||i,25(o iT ; L 29(n)) < C, 

where C depends only on 12 and the final time T. Then 

~{\\u T - il 11 £oo ( 0 ,T ;L 2 )) + l|V ( U T - u) ||l,°°(0,T;L 2 (n)) + l|V(u T — U) \\l 2 (Q,T-,L 2 {0)) ) < C 

for all t £ (—To, To), t / 0, where C does not depend on r. 

Proof. We can rewrite the norm on X = C 1 (0, T; L 2 (Q)) H T/^O, T; 17q (12)) (cf. Proposion 
12 .21) as 



*e{+,-} 


where 


Mix; •“ ( II Hi 00 (0,T ;L 2 (£li)) + ll^^h II Z/°°(0,T;L 2 (n;)) + ll^“* llL 2 (0,T;L 2 (r2;)) 


1/2 


By employing assumption (14.16[) . we can modify estimate (14.14|) by changing the second to 
last line as follows: 


(4.17) 


I (fiiVi)x* x \ 

l ’ L 

< =1 It — l|l,°°(Q i )(l + «o)||flj||L 2 (0,T;L 2 (f2 i ))ll^ T |l^ 2 (0,T;L 2 (n i )) 

+ (l^r — l|L°°(Qi)^l + l^ r _ ^lL°°(Oi)(l + Pl))(-\\^ u i\\L 2 ( 0 ,T-,L 2 (ni)) 

Q 


+ &(1 — ^)l|Vli[|| L 2( 0iT; L2(Q.)))||Vl)[|| i 2( 0]T;i 2(Q.)) 

+ bd(q\A r — 7|£oo(q.)(1 + | A t — I\loo m ) + I A r — 

+ I^T — l|L°°(ni)/3l )l|Vtlj lli2 9 ( 0 iT; L29(Q.))||Vl)j ||L 2 (0,T;L 2 (fi;)) 


+ | 7 t — 1 | 


2k 


T°°(ni) ^ ll^r II L 2 (0,T -.H 1 (12;)) ll^rill, 00 (0,T;L 2 (f2;)) • 


This further implies that 


l*< 


*€{+,-} 


lift 


’ V i)x*Xi 


< C(\I T — 1|l°°( Q .) + \A t — I\ La o^.))\\vJ ||^, 
where C > 0 does not depend on r, from which we can conclude that 

IKIIx — C(l It — l|L°°(ni) + I A t — 7|loo(q.)). 
The assertion then follows by applying Lemma 14. II 


□ 
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5. Auxiliary results 


In order to calculate the shape derivative of our cost functional, we will need to employ 
the two forthcoming propositions. The assertions correspond to hypotheses (H4) and (H3) 

in [Ti] , 

Note that, since 

^ ll^^ T II(0,T;L°°(f2)) = ll^rVK o F r ) ||z,°°(o,T;L°°(0)) 

< (1 + To|T)/l|x,oo( f2 ))||Vu T ||l,°o( 0iT; x, 0 o(n)), 

if we assume that 

{W 2 ) ||Vii||x,oo( 0 )T ;i<x>(n)) < C, 

where C depends only on the (fixed) domain 17 and final time T, we also know that 
l|V« T || L°°(p,T;L°° (o)) is uniformly bounded for |r| < To- Then also condition (j4.16|) holds. 

Proposition 5.1. Assume that the coefficients in the state equation satisfy and q > 2. 
Let hypotheses {TL\) and {Tif) hold. Then 

Jim ~{{E(u t ,t) - E(u,t )) - (E(u T ,n + ) - E(u,Pt + )),p)^ % = 0 

holds for the adjoint state p. 


Proof. We begin by calculating the difference 

\e(u t ,t) - E(u,t) - ( E(u T ,n ) - E(u,n)),p)x^x 

=\ Y, (ii+ih+nh+iv^ 


where we use the following notation 



- uf) — 2 ki(uj — Ui)uij (I T — 1 )pi dx ds 
1)(—2 kiiifpi + (1 — 2kiuJ)pi) dxds 



1 )pi dx ds, 


if = 




A r - /)(—V(tt[ - uf) + bi{ 1 - 5i)V(uf - Ui )) 
Qi 

1 


+ (I T - 1)(—A r V« - + bi{ 1 - 5i)A r X7(ul - u i )))-A T 'Vp i 

Qi 


+ - Ui) + bf 1 - 5i)V(ui - u^fAr - /)Vpj| dxds, 

ni i = jf jT Mi{(|A r Vi£|«- 1 A r Vi£ - \A T \7ui\ q - l A T S7u?j • A T V Pi I T 

~ (|VuJ] ,_1 VuJ' - |Vui| ?_1 V«i) • Vk} dxds, 
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f T f 2ki 

IVi = - / -Ui)(uJ+ Ui)(I T -l)pidxds, 

Jo Jfli 


i € {+,—}. Thanks to hypothesis (Hi), we can estimate the first integral as 

\Ii\ < {^11^1 - ^i||L 2 (0,r;L 2 (n))(2fe||'wI||L 2 (0,r;L° o (0 i )) I |Pi 11 L°° (0,T;L 2 (fl)) + (1 + a o) Wfii IIl 2 (0,T;L 2 (T2)) ) 


2k.. T 

+ t ik 


~ u i||L 2 (0,T;L 2 (a))ll^i||L° o (0,T;L° o (Q i ))l|Pi||L 2 (0,T;L 2 (n))}l^r ~ 


with ao defined as in (12.71) . Integrals //* and /V* can be estimated in a similar manner 
and, by employing the uniform boundedness of A r , the Holder continuity result given in 
Proposition 14.21 and properties of the mapping F T , we can conclude that 

— 'y ) + JJj + IV?) > 0, as t > 0. 

T 

In order to show convergence of the remaining terms to zero, we first rewrite IIli as 

r T 

I Hi = 


(5.2) 


Mi{(|A r V<| 9 - 1 A r V<-|A r V«i| 9 - 1 A r Vu i 

- (IVuJI^VuJ - IVu^-Vu*)) • A r V Pl I r 
+ (|Vu[| 9_1 Vu[ - |Vui| 9_1 Vui) • (I t A t - 7)Vpi} dxds. 


By employing inequality (12.111) (with r\ = 0) and hypothesis (H 2 ) we can estimate the last 
line as follows 

[ [ Mi(|Vu[| 9 _ 1 Vft[ - |V«i| 9 - 1 Vui) • (I t A t - I)X7 Pi dxds 

Jo JfLi 

< C q b5 f T [ \V(u^ -UiMVUill- 1 + \VuU q - 1 )\(I T A T - I)V Pi \dxds, 

Jo Jtli 

— Cq 6 d||V(ftj — ili) llL 2 (0,T;L 2 (n)) (ll^^h Hl oo (0,T;L oo (£2)) 

T ||Vu, 11 £,00 (o,T;L°°(n))) — II lli 2 (0,T;L 2 (fi ))j 

which, after division by r, tends to 0 as r -> 0, due to Lemma 14.11 Proposition 14.21 and 
uniform boundedness of ||Vu 7 ||/ / oo(o i 'r;L° o (r 2 ))- If remains to estimate the first two lines in 
() 5.2 1) . We will first rewrite them using the representation formula (|2.8j) as 

[ [ b i 8i(\A T VuJ\ q - 1 A T VuJ - \A T S7ui\ q - l A T \/ui 

Jo Jfli ' 


ArVpJr 


- (|Vu[| 9_1 Vu[ - |Vui| 9 _ 1 Vu*)) • A T X7piI T dx ds 

[ [ biSJ(A T - /)V(u[ - iii) [ \A T Viii + aA T V(u T i - iii)\ q ~ l da 
Jo Jtti 1 Jo 

+ (V(iiJ - iii) J [\ArViii + aA T V(u\ - ti *)! 9 " 1 - | Vu; + aV(uJ - iii)\ q ~ l ^)da 
+ (q- 1) f (c(A T \7ui + aA r V(iiJ - iii),A T V(uJ - in)) 
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— £(V«i + aV(iiJ — iii), V(uJ — Ui))^ da) • A-Vpj/ r j dxds. 

with C as in (12.91) . For the first line on the right hand side (divided by r) we immediately 
have convergence to 0, thanks to Lemma 14.11 uniform boundedness of A r and Proposition 
14.21 For the remaining terms we further have, due to inequalities (12.111) and (|2.10l) 

f f V (uj — iii) f (\A r Vui +aA T V(uJ - 
Jo Jvti Jo v 

— |Vhj + aV(u[ — in) I 9 ” 1 ^jda ■ A T WpiI T dxds 
<C q \\\7(ul - hi)|| L 2 ( 0]T;i 2 ( Q .))|A T — I\L°°(n i )(\\'V'u-i\\ L °°(o,T-,L°°(n i )) 
d" ll^C^i ~~ 'di) lll/ oo (0,T;L oo (Q i ))) ((^ d~ @1 ) 11^^* II L°° (0,T;L°° (f2i)) 

d- (1 + Pi 2 )l|V(u[ - '“i)llz / oo 2 ( 0r; Loo ( - n .))) ai/dl||Vpi]| L 2(o iT; L2(Q.)). 
Furthermore, by employing estimate (12.141) with 77 = 0 we obtain 

(<?-l) [ (C(A T Vui +aA r V(ui - Ui),A T V(ul - iii)) 

Jo 

— C(Vui + crV(h[ - Ui),V(uJ — iii))) da 
< C q \A r - I\ L oo {Qi) \V(uJ - «i)|{(|V«i| + |V(u[ - iH)\m- 3 + l)(|V«i| 9 - 3 
+ |V« - Ui)\ q - 3 )P! + (|Vu,r 2 + |V« - h,)| 9 - 2 )(l + ^i)(|V*i| 

+ |V(u[ - «i)|)}, 

valid for q > 2, from which, by making use of Proposition 14.21 Lemma 14.11 and hypothesis 
(7ii), we finally have \llli —» 0 as t —>• 0 and therefore Proposition 15.11 holds. □ 

For the second property to hold we have to assume that p is slightly more than W 1,2 regular 
on the subdomains, i.e.: 


(Ti 3 ) P\Qi G L 2+£ (0,T; VF 1 ’ 2 ^^)) for some e > 0. 

Proposition 5.2. Let q > 2 and assumptions (Em on the coefficients hold. Assume that 
hypotheses (TLi)-(TL 3 ) are valid. Then 

lim —(E(u T , n + ) — E(u, fi + ) — E u (u, Ll + )(u T — u),p) y* y = 0, 

r—>■ 0 T A 

where p is the adjoint state. 


Proof. We have 


(5.3) 


(E(u , n_j_) E(u, 12+) E u (u,Ll^-)(u u) , p)j> 



«i) + « 



+ E f T f^ifiVuir^ul-lVui^Vui-lVuil^Viul-u^ 

*={+,-} ■ /o Jn V 
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— (q — 1)| \7ui\ q 3 (Vtq • — Ui))\7ui'j ■ S7pi dx ds. 

The first sum on the right hand side can be estimated as follows 


E 

*£{+,-} 

! E 

*€{+, — } 


r T r 2 k t 

( u i ~ Ui)(uJ - Ui ) + (u[ - iiiY ) pidxds 


10 JO; 


2 h 


A,; 


A* 


« - Ui)(ul - Ui)pidx 


l T 


J 0 


+ 


E fj a ft 

.Jrr i jo Jo 


— Ui)(uJ — iii)pi dx ds 


*£{+,-} 


2k 


^ ^ -r i ( C '2i,L 4 ) 2 |l U <’ ~~ II L°°(0,T 1 (O;)) ll^r “ '“illL° o (0,T;L 2 (ni))ll^*lli 2 (0,T;_H' 1 (ni))’ 


ie{+,-} At 

since (u — u T )\t=o = (it — u T )|i=o = 0 and p\t=T = P\t=T = 0. This expression, upon division 
by r, tends to 0 as r —>■ 0, due to Proposition 14.21 The second sum in (15.31) can be rewritten 
with the help of formula (12.81) as given below 


ie{+,-K° 

— (q — l)|Vui|' ? ” 3 (Vui • V(ti[ — u)i)Vu^j ■ Vpidxds 

f f biSi f (|Vuj + crV(u[ — fiJI 9-1 — iVuil 9-1 ) da V(uJ — «*) • Vpj dx ds 

r_l_| J 0 J fli J 0 

+ ^ j j bi5i(q-l)J (c(Vui + aV(u T - in), V(uJ - iti)) 


ie{+,-} 

— £(Vuj, V( 1 J — iti))} ■ X7pi da dx ds 
:=/ + //, 


where £ can be estimated as in (|2.14f) . By employing inequality (12.111) we obtain 

|/|< V [ T [ b6C q \V(uJ-iti)] 2 - 71 [\\Vui + aV(uJ-u^l + lViiiiy-^dalVpildxds. 
^r_| | Jo J£li Jo 


Making use of Holder’s inequality and hypothesis ('Hi) results in 


\i\<bsc q Y T [ l y ( 

Jo Jo 


itj - Ui)\ 2 -^\Vuir 2 ^ + IV(u[ - dxds 


*e{+ 

<bSC q Y II'^ 

*£{+,-} 

+ iiv«-ujiir4+v 


U i ^)Hl 2 (0,T;L 2 (O)) 


CII II 9-2+77 

111 V IIL°°(0,T;Z,°°(0;)) 


L~(0,T;L«(O))) II'H L § (0iT;L | (n)) - 


Here we can choose rj = ^- 7 . Recall that A||V(tJ —uJll^o^LjO)) is uniformly bounded for 
|r| < ro, r J 0, due to Proposition 14.31 and that thanks to hypothesis (H 2 ) we have uniform 
boundedness of || Vu1||j,cx ! (o,T;I/ x> ( 0)) as well. This means that we can achieve that / upon 
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division by r tends to zero as r —> 0. 

By employing inequality (12.141) . we get the estimate for II: 

\II\<bSC q V [ T [ |V(u[ - ix)| 2 ^{|Vi*i| 9_3+TJ + |V« - Ui)\ q ~ 3 +r, )(\Vui\ + |Vu[|) 

+ |Vuj|' ?-2 (|Vhj| T? + |V-^XT)} 

<bSC q ^ l|V(u — Wi)|| i 2^ 0jT;i 2(Q.))|(||Vui]|| (OO ^; ioO ( n .)) 

* 6 {+,—} 

+ II “ «Ollloc 3 ( o"T;L-(a i )))(ll V ^ll^(0,T;L-(a i) ) + II Hl°°( 0,T;L°° (n ;))) 

+ ll^ 7 ^llloo( 0 jT;L oo(Q))(||Vu||^ oo ^ 0iT . ioo ^.^ + ||Vhj ||^( 0 , T ; L .(, 0 ))}l|V K |l L a (or . L | (ni )), 

with r) = Upon division by r, due to Propositions 14.21 and 14.31 the right hand side tends 
to zero as t — ^ 0. □ 


6. Computation of the shape derivative 


Let u T , u satisfy E(u t ,t) = 0 and E(u,Q + ) = 0, for |r| < r 0 , r € R. In that case 
u T = u T o F r is the solution of E(u r , fi+ iT ) = 0. We then have 

rT 


dJ(u,Q+)h = lim — [ [ (j(u r )I T — j(u)) dx ds 

T ^° T Jo Jn 


[ [ (j'( u ) liiri U _ U + j(u)dWh) dx ds , 

Jo Jn 


t— >■ o r 

where we have used (similarly to Lemma 2.1, p3]) that 

1 r T 


lim 

t —>0 T 


< lim — 

r—>0 r 


I /(*(«') - jN - f(u)(u T - u)^jI T dxds 
f f (u T — u) 2 I T dx ds 

Jo Jn 


which follows from Proposition l4.2l and the fact that / r is uniformly bounded for r £ (—To, To). 
By employing the adjoint problem (12.161) and then proceeding as in the proof of Theorem 
2.1, [14], we obtain 

/ / j'(u)(u T - u)dxds = (E u (u, Q. + )(u T - u),p)^ % 

Jo Jn 

= - (E(u T ,Q + ) - E(u,Q+) - E u (u,Q + )(u T - u),p)x*,x 
~ (E(u t ,t) —E(u, r) - ( E(u T , fi+) - E(u, Q + )),p)^ i ,^ 

=o 

- ( E(u,t ) - ^(w,0),p)^* ^. 

=o 

The second and third line divided by r tend to zero, as r —> 0, on account of Propositions 
o and 15.21 and we are left with 
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= - lim ~(E(u,t) - E{u,ti),p)x* f 


t->o r 


This limit, representing the (artificial) time derivative of (E(u, is typically com¬ 

puted by transforming the expressions E(u , r) and p back to E(u o F~ x , fI +)T ) and p o F^ 1 , 
and then making use of differentiation rules for u o F^ 1 and p o F (T 1 (see Examples 1-5 and 
Lemma 2.4, M)- However, these rules assume H 2 differentiability in space of the primal 
and the adjoint state, which is too high of a requirement in our case. Instead, we continue 
with calculating the difference 

(E(u,t) - E{u,ti),p)x*'X = I + + I -i 
where the two terms on the right hand side are given by 


Ii = 


f ( {^-(1 - 2kiUi)uiPi{I T - 1 ) 

J 0 Jfli 1 A i 


+ — ((A r - I)Vui ■ X7pi + A r Vui ■ (A t - I)Vpi + A T Vui ■ A T X7pi(I T - 1)) 

Qi 

+ 6,(1 - 8i)((A T - I)Vui ■ Vpi + A r Vui ■ (A t - I)Vpi + A r X7ui ■ A T X7pi(I T - 1)) 
(6.1) + bi5i((A T — I)Vui J \Viii + a(A T — I)'Vui\ q ~ 1 da ■ A T VpiI T 

+ (q — 1) [ |Vu, + a(A T — I)Vu,| 9_3 (Vu, + a(A r — I)Vu,) 

Jo 

■ (A t — I)Vui(Vui + a{A T — 7)Vtt,) da ■ A T \7piI T 
+ I VUiP-'VUi • ((A t - I)VpiI T + (J T - l)Vpi)) - ^{ui) 2 Pi{Ir - 1)} 


dx ds, 


i G {+, —}, and we have employed the formula (12.81) to represent the difference \A r 'Vui\ q 1 A r 'Vui~ 
| Vu,| 9 _ 1 Vuj. Dividing (16.11) by r, passing to the limit and utilizing Lemma 14. 1 1 yields 

lim - (E{u, r) - E(u, 0),p)x* j 

T 

= V / / {-(-Vuf + b i {l-8 i )Vii[ + b l 5 i \Vu i \ q - l VuJ)(Dh T Vp i + DhVp i ) 

(6.2) + bi8i(q — l)|V'u,| 9_3 (Vuj • (— Dh) T \7iii)(Vui ■ Vp,;)| dx ds 




ie{+,-} 


- 2 kiUi)uiPi + —Vui ■ Vpi + 6,(1 - 8i)\7ui ■ Vpi 
Qi 

2h 


+ bi8i\Vui\ q 1 Vu, • Vpi - r-{ui) 2 pi f div/i dx ds, 

and we can now express the Eulerian derivative: 

Theorem 6.1. (Weak shape derivative) Let q > 2, uq,u\ G W 0 1 , 9 + 1 (H) ; and assumptions 
(13.11) on coefficients hold. Assume that {Tli)-{Tlz) are valid. Then the shape derivative of J 
at fl + with respect to h G C' 1,1 (D, M d ) can be represented as 

dJ(u,n + )h= f [ \ (-Vu T + 6(1 — 5)Vu t + b5\Vu\ q - 1 'Vu T )(Dh T 'Vp + DKVp) 

Jo Jn i Q 
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(6.3) 


+ bd(q - l)|Vu| 9_3 (Vu • Dh T \7u)(yu • Vp)} dxds 

— / {--(1 — 2 ku)up -\—Vu • Vp + 6(1 — 5)Vu ■ Vp 

io isi 1 -'' Q 

2k I 

+ bd\\/u\ q ~ 1 Vu ■ Vp — —{u) 2 p — j(u)j di vhdxds. 


Note that the integrals in (16.31) are well-defined thanks to hypothesis (Hi), and for them 
to be well-defined hypothesis (Hf) is actually not necessary. 

Theorem 16.11 gives us the shape derivative of the cost functional in terms of the volume 
integrals, which is in p] regarded as a weak shape derivative. However, an obvious advantage 
of the volume expression of the shape derivative is that it allows for a lower regularity of 
shapes as well as the lower regularity of the primal and the adjoint state. Recently there have 
been suggestions that the domain representation is also advantageous in terms of easiness 
of computation and numerical implementations (see, for example, EH], [22]), especially in 
the case of transmission problems where shape derivatives given in terms of the boundary 
integrals contain jumps of functions over the interfaces, which is numerically a delicate task 
to perform. 


6.1. Strong shape derivative. In order to express the shape derivative in the form required 
by the Delfour-Hadamard-Zolesio structure theorem we would have to apply Green’s theorem 
to the last two lines in (16.21) . which is not allowed since u and p are not sufficiently regular. 
However, it turns out that if the domains are sufficiently smooth and Vu is bounded in 
L°°(0, T;L°°(H)), the state variable exhibits H 2 regularity on each of the subdomains. This 
result together with an assumption regarding the regularity of the trace of Vu| n and Vp| n 
on T makes expressing the shape derivative of the cost functional in terms of the boundary 
integrals possible. Let dQ and T = be C 1,1 regular. We utilize the following result (cf. 
Theorem 2, my- 

Theorem 6.2. Assume that q > 1, q > d-1, uok € H 2 (fU «o,«i G W 0 ll9+1 (fi), 
and that <9f i and T = are C 1,1 regular. Let u be the weak solution of (|2.3|) . If u € 

l / L 1,oo (0,T; W 1,0O (fi)) and \\\/u\\ L oo^ 0T . Lao ^^ is sufficiently small, thenui € iL 1 (0,T; iL 2 (H,)) ; 
i € {+, —}• 


To be able to express the shape derivative in terms of the boundary integrals over T, we 
first impose additional regularity hypotheses on u and p: 

(Hi) tr£‘Vp € L 1 (0, T; L 1 (T)), i € {+,-}, 

(fts) tr^Vu € L°°(0, T; L°°(T)), * € {+, 

where ti'p ! Vu and tr^Vp stand for the trace of Vu|n i and Vpk, respectively, on T. Hy¬ 
potheses ( 7 ^ 4 ) and (fHf) will ensure that the forthcoming boundary integrals are well-defined. 
Note that they do not follow from the previous hypotheses and regularity results, 

partially due to the lack of an appropriate trace theorem in the limiting L°° case. 

Next, we introduce sufficiently smooth in space approximations of the adjoint state in 
i7 1 (0, T; i7 1 (Hj)). Fix i € {+,—}■ Let {pi, m }m =1 c iL 1 ( 0, T; C°°(fli)) be a sequence that 
converges to pi in i7 1 (0, T; iL 1 (Hj)) and such that pi tTn = Pi on (cf. Theorem 3.42, [6]). 
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We can then approximate (16.2D as 


I™ \ i E (u, t) - E(u,0),p)x 

E . 

*£{+,—} 


x*,x 


1 


-(-Vuf + 6,(1 - <5,)Vhf + 6 i ^|V?x i |' ? - 1 Vuf)( J D6 T Vp i 

'O do, 1 ft 


(6.4) + DhVpi^m) + bi5i(q - l)|Vh,| 9 3 (Vh, • (-Dh) T Vui)(X7ui ■ Vp,, m ) j dxds 

pT C 

4“ ^ ' / / "I t - ( 1 m H Vnj • m -I- 6,(1 d,)Viij • Vp, 5 

iG ^_ } 4o ./n. ( >• Ai ft 

2jo. .. 

+ 6j^|Vrt,| 9_1 Vtt, • Vp,, m - ^(hi) 2 ^,™ f div/idxds + Ri(pi,Pi, m ), 

Ai > 


where the error term is given by 

r T 


1 


-(-Vuf + bi{l-6i)Vuf 
ft 


Rl(Pii Pi,m) — ^ ' / / { (" 

,£{+,-} ^ * 

+ 6j<5,|V«j|' ?_1 Vt(f)(D6 T V(p, -Pi, m ) + DhV{pi - pi !m )) 

+ 6A(g - 1)| Vhj| 9 “ 3 (Vu, • (-Dh) T Vui)(yui • V(pi -ft,™))} dxds 
-T /• ^ ^ 

^ ^ / / ] T - (1 2kiUi)Ui(pi Pi rn) H ^(jPi Pi,m) 

iejE)-' 0 • /n ‘ lAi e ‘ 

H“ ^z(l ^()Pz Pi,m) H“ ^Q?z Pi,m) 

- (h,) 2 (ft - Pi t m) } div/i dx ds. 

Since n, and ft jm are sufficiently smooth, we are allowed to employ Green’s theorem in (16.411 . 
This will cause the terms containing Dh (not included in i?i) to cancel out, and we arrive at 


iim - (E(u,r) - E(u,0),p) 

r-fO T 


X*,X 


E 


*e{+,-} 


/O JdQi 


/ (1 2kiUi)UiPi^ m ~\~ V Ui Vp^,m 

2k- 'i 

+ 6,(1 - dj)Vhj • Vft,™ + 6,d,|Vtij| 9_1 Vui • Vft,™ - -^(-iij) 2 ft,™ j h T rndxds 

pT P \ 

- Y1 / / {tt( 1 “ 2k iUi)ui{\7pl m h) + —Vui • V(Vp^ m h) 

,e{+,-K° ^ ^ 

+ 6,(1 - dj) v ^i ■'V(^ / Pi,m h ) + 6,dj|Vh,| 9_1 Vhj • V(Vp^ m h)| dxds 

pT P ~\_ 1 . 

~y2 {W—(1 — 2 kiUi)pm t i(VuJh)-\ -Vft,™ • V(Vufh) 

is{+,-} Jo 4n, L A, ft 

- 6,(1 - di)Vp m ,i • V(Vrtfh) - 6A (G Mi (Vft,™))' • V(Vufh)} dxds + Ri(pi,Pi, m ). 
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This expression can be rewritten as 

lim — (E(u, t) — E(u , 0 ),p) 


t— >0 r 


A'*, A' 


E / / |t-( 1-2 kiUi)uj,pi 

r~ , io Jaiii 


*€{+,-} 


H-Vuj • Vpi 

Qi 


2 k- ^ 

+ bi{ 1 - • Vpi + bi5i\S7iii\ q ~ l S7iii ■ Vpt - -^(ikfpi j h T mdxds 


£ 

*G{+,-} 


n 
/ 0 


T-(l “ 2 kiUi)ui(Vpl m h) + —Vui ■ V(Vpl m h) 
A i Qi 


+ /)i(l - 8 i)Vui ■ V(Vpf m h) + M.IVitil 9 1 Vu; - V(Vp^ m /i)| dxds 


£ 

*€{ + , — } 


/o 


{T-(l - 2 kiUi)pi(yujh) + —Vpi • V(Vuf/i) 

*• A,; Pi 


- &i(l - <5*)Vpi • V(Vufh) - 6A (G^Vpi))- • V(V«f/i)} dxds 
T R\(piiPi t m) T -R 2 (PiiPi,m) T ^3 (Pi : Pi,m) i 

(the first and the third sum are written in terms of pi plus the error R 2 + R 3 ) where the 
approximation error terms are given by 

rT r , 1 

- 2 kiUi){pi - Pm,i)(yujh) 


R 2 (Pi,Pi,m)= V / / (v^ 1 

1 _\ 0 J Ui i 


*e{+,—}' 

+ — V(pi -pm,i) ■ V(Vttfh) - 6j(l - 5i)V(pi -p m ,i) - y{yujh) 

Qi 

- bA (G Ui (V (jpi -Pi,m))) m • V(Vitfh) j dxds, 

RsiPiiPi^m) = ^ ^ ll\ 2kiUi)Ui(j)i Pi,m) H ^ (Pi Pi,m) 

iG | + _j JO JdQ* 1 

H“ ^z(l ^z)^^z ' ^ (p* Pi,m) H - Mz | V* V (Pi Pi,m) 

- ^(hi) 2 ( P i -Pi,m)| /i r Ui dxds. 

Remark 6.1. If u £ Vk 1,oo (0, T; TT 1,00 (f2)) ; we are allowed to use <f £ L 2 (0, T; Rg (0)) as a 
test function in (11.31) . To see this, let 4>j £ L 2 (0, T; Co°(fZ)), 4>j -A cf in L 2 (0,T;R o 1 (fi)) as 
j —> 00 . Then 

rT 


f [ { — (1 — 2 ku)ucf)-\ —V-u • V0 + 6(1 — <5)Vti • V0 
•7o io A Q 


2 k 
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< 


A 


((1 + a 0 ) ||w||l 2 ( 0 ,T;L 2 (Q)) + 2fe(C^i L 4 ) V^T'|| , u|| ^oo ( f2 ))) II ^ *A? ll-I/ 2 (0,T;L 2 (f2)) 


,1 


+ ( - ll^ u llL 2 (0,T;Z, 2 (n)) + b(l ~ ^)l|V'u|| i 2 ( 'o ) T;L 2 (n)) 

— 

+ M||Vu||| /OO 1 ^ O)T . Loo ^^||V'u|| Z/ 2 ( O)T;i 2 (Q)))||V (0 — (pj)\\L 2 ( 0 ,T-,L 2 (Cl)) 

—>• 0, as j —>• oo. 

Note that X7ufh € L 2 (0, T; iY" 1 ^)) and Vpf m h € L 2 (0, T; however functions 

I V?l.T 11" f'l/jY- tA T £= O , 

<KM):= 


Vrt^(x, t)/i(x) i£Sl + 

V/(j;,i)li(i) x € A2_ 


> r .X f ^Pm,+(x,t)h(x) xeQ + 
\ Vpm_(x,t)h(x) x € 


do not have to be continuous across the boundary T and we cannot use them directly as test 
functions in the weak formulations of the state and the adjoint problem. We can instead 
employ the two-donrain weak formulations which results in 

lirn ~{E(u,t) - E(u,0),p )** % 

M 0 T ’ 


E 

*€{+,-} 


I!' 

1 0 Jdfli 


j t (1 - 2 kiUi)uiPi + —ViLi ■ Vpi 

^ Ai Qi 


2 ki 


+ bi{ 1 - <5i)Vwi • Vpj + Mi|Vai| 9 1 Vft i • Vpi - ^(wi) 2 ? 3 *} 


/i T nj dx ds 


E 

*£{+-} 


f T J 


{ — W" 1 +bi(l-Si)^ +bi5i\Vui\ q 1 ?^-\ (V pj m h ) 
t Pi orii oni am j 


dii. 

* <9n 


_i du t 


dx ds 


V r [ \-^-b i (l-6 i )p--biS i (G Ui (Vp i )- ni y}(VuTh)dxd8 
fr' x Jo JdQi 1 ft. dn* 5ni J 


*g{+,-}' 


y2 f(Ui)(Vufh)dxds + Ri(pi,Pi, m )+R 2 (pi,Pi,m) + R3(Pi,Pi,m)- 

v , Jo Jfti 


*G{+, — } 

Finally, this can be rewritten as 

lini — (E(u,r) - E{u,Q)),p)x* x 

r-iO T ’ 


pT C \ 

Tl t-(1 - 2 kiUi)uiPi H-Vrij • Vpi 

r , ! J 0 YaOi *- A i Qi 


*€{ + , — } 


2k- 'i 

+ 6j(l - 5i)Vui ■ Vpi + bi5i\S7Ui\ q ~ l S7in ■ X7pi - ^-{urfpi } h T mdxds 

Ai J 


E 

ie{+,-} 


n 

/ 0 Jdfli 


E / / 

ie{+,-} J° J 9Qi 


+ M 1 “ + bi5i\\7ui\ q (Vpf/i) deeds 

l Pi OUi OHi OUi > 

{ — 7^ - &i(l - ~ Mi (G Ui (Vpi) ■ UiY | ( Vufh)dxds 
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yZ [ [ j'(ui)(Vufh)dxds 


-|- R\(pi, Pi t m) T R-2 (Pii Pi,m) T Rs^Pii Pi,m) “I - R& {Pii Pi,m)i 
(the second sum is written in terms of pi plus the error R 4 ) with 

1 dui diii 


R^ijPii Pi,m) — /* ' / / 

i£ { + _}J0 JdQi 


f)qt 

+ 6,(1 - 6,)^ + MilVAil^V 

Qi orii orii orii 


l - f v(p, - Pi : m) T h dx ds. 

i s 


Next, we want to show that 

R{jPiiPi,rn) •— RlijPiiPi,m) T R2(PiiPi,m) T RsiPii Pi,m) “1“ R&{]?ii Pi,m) 
tends to zero as m —> 00 . We will focus here on the estimates for the boundary integrals. 

Since pi — pi >m = 0 on T, we know that Vr(p, — Pi,m) = 0, where Vr denotes the tangential 
gradient. This further implies that 


(6.5) 

( 6 . 6 ) 


\71 \| t 9(pi Pi m) /, 

V(p, -Pi,m)\r ■ h = - — - (h ■ Hi), 

d(pi - Pi,m ) dui 


V(Pi~Pi,m) |r • Vuj| r = 


drii drii 

Due to (16.51) and the fact that h = 0 on dQ- \ T, we can estimate 

diij 
dru 


Ri(Pi,Pi,m) ~ J2 J 0 J m L i d Sni +bi{l 


*£{+»-} 

+ b i 6 i \Wu i \ q - 1 ^] 

<ih| L cx)(r) y 
*e{+,—} 

19-1 


-‘^1 a{V '- pi ’ m) .h T n,d x is 


dm 


(I 

dm 


dm 


I/ 2 (0,T;H 1 / 2 (r)) 


+ Wi-« 8i “ 


dm 


L 2 (0,T;_H 1 / 2 (r)) 


du, 


+ 6i<Ji||V«i|| Lao(0iT;ioo(r)) || ^ 


) 

<9(p, - Pi,m) 

L 2 (0,T;H 1 /2(r))/ 

dm 


<c\h\ L oo( T ) y 

*€{+, — } 


fl 



5rij 


■ 0, as m —>• 00 . 


L 2 (0,T;if 1 / 2 (r)) 
dm 


+ 6,(1 — Si) 


din 


L 2 (0,T;_H _1 / 2 (r)) 
L 2 (0,T;iL 1 /2(r))y 


L 2 (0,T;_H 1 / 2 (r)) 


dm 

II Pi ~ PhmW^iO^H 1 ^)) 


dui 


Here we have made use of the fact that since m G i6 1 (0, T; H 2 (Qi)), we have € 

H 1 (0, T; H 1 / 2 (dfli)) (provided that Qi has a C 1,1 boundary, which we have assumed). Sim¬ 
ilarly, by employing (16.61) . we obtain 

rT 

R3(pi,Pi,m) = ^ ' 


iG{+,-} 


f 1 (, n , v ( \ , 1 dui d(pi-pi t r 

I I | \ (1 ^kiUijUiyPi Pi,m) 

/0 L A i 


din d(pi - pi, m ) | 

T 6,(1 Oj) rj o + biSi Vuj 

cm, cm. 


Qi dm dm 

—X dlli d(pi Pi,m) 

dm dm 






































30 


V. NIKOLIC AND B. KALTENBACHER 


- ^(iii) 2 (pi - Pi,m) \ h 1 riidxds 


< 


ie{+,-} 


(L 

dm 

\Qi 

dm 


L 2 (o,T-,m/ 2 {r)) 
19-1 


+ h{ 1 — Si) 


dm. 


l^li°°(r) ^2 {^ _ l|iii||L 2 ( 0 ,T;L 2 (r))l|l — 2fc;/Uj|| I ,oo(o j T;L“(r))l|Pi ~ Pi,m\\L 2 (0,T;L 2 (?)) 

dm 

L 2 ( 0 ,T;_ff 1 / 2 (r)) 

II Pi ~ 

+ II III/ 00 (0,T;L°°(r)) IIlij IIL 2 (0,T;L 2 (r)) ||Pi — Pi,m||.L 2 (0,T;L 2 (r))} 

• 0, as mm- oo. 


,1 t lly 11 g— 1 II duj 

+ OiOiW VRi|| Z/ oo( 0iT;i oo( r ))|| — 


L 2 ( 0 ,T;_ff 1 / 2 (r)) 


2h 


Altogether, this means that 


1 


lim - (E(u, t) - E(u,0),p)^ * 
T-* o r ’ 


^ , Jo idfli ^ ^ 


— 2k, 


*£{+>—} 


H-Vrtj • Vpi 

Qi 


2k- 'i 

+ 6i(l - Si)Viii • Vpi + M*|Vzii| 9_1 Viit • \7pi - ——{ui) 2 Pi \ h T riidxds 

A; ) 


+ bi{1 - *>£ + W|V4i|,_1 ^} ^h)dxds 


{“T 5 - ~ M 1 " - ''.A (G ul (V Pi ) ■ Hi )'} (Vuf h)dxds 

^ 4/4 c j iii (jiii / 


V 

;eh7-} J° J dfli 

^2 f f j'( u i)(y u Th) dx ds. 


r\ L\ fc / <~v 

Qi dm dm 


*€{+,-} 


/o JQi 


Assume that Ud € L 2 (0, T; We can utilize the Stokes theorem and the fact that 

h = 0 on d£l- \ T to acquire 

/2 / / d( u i) div/idx ds + / / j'(ui)(S7uJh) dx ds 

s J \ ^ 0 J zE{H } ^^ ^ 


*€{+,-} 




/in* dx ds, 


*€{+,-}• 

which leads to the shape derivative given in terms of the boundary integrals 


dJ(u 


, 0_|_)d = — f [ —(1 — 2 ku)up-\ —V« • Vp + 6(1 — dlV-d • Vp 
do dr 11 A £» 

2k n 

+ 6(5|V'd| ,? “ 1 Vu • Vp — {u) 2 p h T n + dxds 

A 

+ [ [ f(-Vu + 6(1 - <5)Vh + 6d|Vh|' 3_1 Vh) • n+(Vp • h) 

Jo dr U- Q 


dx ds 
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+ £ f ~ 6(1 ~ w - b ^ G ^ v p^y) • n +( vu • h ) 


dx ds. 


Here we have used that [[(it — Ud) 2 } = 0. The expression for the shape derivative can be 
slightly simplified. For the second and third integral on the right hand side, by employing 
the fact that (see Example 2, [14] ) 

M = M = 0 =>- [xyj = 0, 
we obtain the following identities 

r (ivu + 6(1 - J)Vit + fe(5|Vu| 9_i Vi[) ■ n + (Vp • ft) 

+ K1 -S)p-^ + 66|V»|«-^^)(ft ■ n + ) 


dn+ d n_i 


dn+ dn .. 


g <9n+ dn 

+ (ivu + 6(1 - J)Vli + 66|Vu| 9_1 Vft( ■ n+Vrp ■ ft 

| i 8 !L ^ + M^ + 6 ^ 

g dn + on .| 


<9ra+ <9ra 


+ 


<9ra+ <9nj 


(h • n + ). 


and similarly 


pT p i 

/ / ^-Vp — 6(1 — d)Vp — b5{G u {S/p))'^j ■ n+(Vu • h) 

r>T 


dx ds 


1 

-Q 

1 du dp i <9?i <9p i Lnr7 ,| j chi dp 


p iL an+ on + 

+ 64(g- l)|Vn|«“ 3 (Vn-Vp) 


chr+ chu 


+ 6d|Vu| 5 


<9n+ 9n_i 


du 

2 1 

<9n+ 



(/i • n + ) dx ds. 


Here we have made use of the fact that [[Vpu] = [[Vpp]] = 0. We finally obtain 

Theorem 6.3. (Strong shape derivative) Let <912 and T = 912+ be C 1,1 , uo\fii € 
uo,u\ € H / 0 1 ’ 9+1 (H), q > 2, and let assumptions (13.11) on the coefficients in the state equa¬ 
tion and hypotheses {TLi)-fHf) hold true. Assume that Ud £ L 2 (0, T; idg(12)). The shape 
derivative of J at 12+ in the direction of a vector field h £ C 1,1 (H,M rf ) is given by 

rT 


dJ(u,Ll + )h.= ( [ -(1 — 2 ku)iip -Vu • Vp 

do dr 11 A g 


(6.7) 


2k 

- b(( 1 - 6) + d|Vu| 9_1 )Vu • Vp + —(n) 2 p 

A 

+ + 26(0-6)+6 |V»| 


p chr+ <9n 
+ 6%- l)|Vn| 9 “ 3 (Vtt- Vp) 


du 

2 1 

dn + 



<9n+ <9n+ 
h T n + dx ds. 


The boundary integrals in (16.71) are well-defined thanks to hypotheses ifH/f) and ( TL§ 


7. Conclusion 

We have computed, through a variational approach, the weak and the strong shape deriv¬ 
ative for the cost functional determining the acoustic pressure of high intensity ultrasound 
when focusing is performed by an acoustic lens. 
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Future research will be directed at developing and implementing a suitable gradient based 
optimization algorithm, as well as considering a physically more involved elastic model for 
the focusing lens and thus a shape optimization problem with an elastic-acoutic coupling as 
the optimization constraint. 
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